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Abstract. We perform an exact computation of the grand partition function of a model of 
confined quarks at arbitrary temperatures and quark chemical potentials. The model is inspired 
by a version of QCD where the perturbative BRST symmetry is broken in the infrared, while 
perturbative QCD is recovered in the ultraviolet. The theory leads, even at tree level, to a 
quark mass function compatible with nonperturbative analyses of lattice simulations and Dyson- 
Schwinger equations. In spite of being defined at tree level, the model produces a non-trivial 
and stable thermodynamic behaviour at arbitrary values of temperature or chemical potential. 
Results for the pressure and the trace anomaly as a function of temperature qualitatively 
resemble those of non-perturbative interactions as observed in lattice simulations. The cold and 
dense thermodynamics is also contains non-trivial features, being unlike a gas of free massive 
particles. 


1. Introduction 

Given the difficulty to address the problem of confinement in strongly interacting systems directly 
from its fundamental theory, Quantum Chromodynamics (QCD), several alternative approaches 
have been proposed. One of the main such approaches is that of effective models of QCD, which 
are quantum field theories that possess one or more fundamental aspects of the original theory 
but being nevertheless easier to have some information extracted. 

Regarding the quark sector of QCD, two quite sucessful models in the description of chiral 
symmetry breaking and its restoration at high temperature are the Linear Sigma Model with 
quarks (LSM) [T], and the Nambu-Jona-Lasinio (NJL) Model [2]. In their original formulations, 
these models do not address the issue of quark confinement. Indeed, in both models quarks 
are effectively treated as on-shell quasiparticles. Two possible directions that can be taken 
towards a simultaneous description of both chiral and confinement dynamics are represented 
either by the coupling of the Polyakov loop to quark degrees of freedom (the so-called PLSM 
and PN JL models) (Sj HI El EJ [71E] , or by considering nonlocal interactions between quarks as a 
result of their nonperturbative coupling to gluons, as in the nonlocal versions of the NJL model 
iSEniiiiiiiTiiis]- Following dlEEllinillTlEH!, we consider a third possibility, which is inspired 
by the (refined) Gribov-Zwanziger effective theory for infrared QCD [19LI2UI ITT] , although not 
equivalent to it. 

^ In collaboration with Leticia F. Palhares and Marcelo S. Guimaraes. 


Although it has been long clear that quarks and gluons are confined to hadrons, a definite 
theoretical criterion for confinement is not yet a settled issue. A sufficient condition for the 
absence of isolated quarks or gluons from asymptotic states (i.e., confinement) is assumption that 
they violate reflection positivity [23]. We therefore take violation of reflection positivity as our 
criterion for confinement, following, e.g., uniiis]. In this work, we explore the thermodynamics 
of a quark model in which confinement is encoded in the positivity violation of the quark 
propagator. 


2. A nontrivial solvable quark model 

As evidence from lattice QCD |23| and Dyson-Schwinger Equations [23] studies shows, the zero- 
temperature quark propagator can be quite well parametrized by a momentum-dependent mass 
function compatible with the functional form 
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where mo is the quark current mass. Indeed, the data of |2j| can be well fitted by Eq. 0 , with 
the values A = 0.196GeV^, m^ = 0.639GeV^, and mo = 0.014MeV [Tn|. It is interesting to 
notice that, with such parameters, the resulting euclidean quark propagator 
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displays violation of reflection positivity, indicating the confinement of quarks m- 

The mass function Q can be obtained from the lowest-level quark propagator of the theory 
given by sum of the QGD lagrangian, 
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with the BRST invariant action 
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between the quark fields d and the auxiliary helds ^ and A. The resulting action 


S = Sqcd + S^x + Sm (6) 

has been shown to be renormalizable m, being equivalent to QCD in the high-energy limit but 
radically changing the infrared sector of the theory (with respect to the perturbative picture). 
We interpret the extra fields and their interaction with the quark fields V' as a local way (in 
the sense of QFT) to take into account an effective dressing of quarks by gluons. The auxiliary 
fields of the quark sector of Q, can be straightforwardly integrated out, providing an effective 
theory for the quarks, whose (nonlocal) action reads, in the quadratic approximation. 
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Notice that, in this model, the quark mass function 0 is directly derived from the effective 
action Q, with A = 2 M 1 M 2 . 

Although for now we only shall investigate the model in the quadratic level, it is perfectly 
possible to consider a loop expansion of the free energy [26]. It is also reasuring that the full 
local model Q is a renormalizable QFT, as shown in US). 


3. The partition function 

At lowest order, the theory defined by Q has a quadratic action. Therefore, its grand partition 
function can be exactly calculated using standard techniques of Finite-Temperature Field Theory 
[26] . Besides having a quadratic action, the theory does not correspond to a free theory, a fact 
that can be clearly seen, e.g., from the nonlocal form of the equivalent effective action Q. 

One can easily introduce temperature by compactifying one euclidean direction (convention¬ 
ally the 4-direction), whereas the introduction of the chemical potential is more subtle. Starting 
from the local action Q, one must first calculate the hamiltonian %. By identifying the quark 
number with the charge associated with the U{1) global symmetry transformation 
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one can use Noether’s theorem to calculate the corresponding conserved current density and 
let AA = / d?xj^ be the quark number operator. The resulting grand partition function 


Z(r, n) = Tr exp 
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can be cast in a functional integral form and straightforwardly calculated |18j . The result is 
convinently split into a /x—independent term plus a (T, /x)—dependent term as 
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where we used the standard sum-integral notation, 
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The p = 0 term can be split in a sum of four terms, two of which corresponding to positive 
pressures of particles with complex conjugate masses, one of which of a particle with real mass, 
and one of which with a negative contribution or the pressure. Each of this term can be calculated 
straightforwardly as the pressure of a free gas |26|. After the usual vacuum energy subtraction 
[26] . one finds 
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for the /X = 0 contribution, where 


log Zq = 2NcNfl3V j {(fi + P 2 + (P3- 2 (^ 0 ) 


(13) 

















is the pure vacuum contribution. The quantities (pi {i = 0,1,2,3) are known functions of the 
internal momentum [TH] . 

The /r—dependent term is given by 
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A closed expression for (14) can be found at the zero-temperature limit. Using Cauchy theorem, 
one may write 
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The expressions above allowed us to evaluate the partition function (10) exactly in order to 
compute several thermodynamical quantities. We present our results in the next section. 


4. Results 

From the partition function, one can derive any equilibrium thermodynamical quantities. Let 
us first show our results for the pressure 


P{T,p) = ^logZ{T,p). (18) 

In Figure we show our results for zero chemical potential, with a comparison between our 
model, a massless bag model and a gas of free massive particles with mass Notice that, 

in the low temperature regime, the model result is quite close to that of a free, massive gas of 
mass Mtfir = 0.467 GeV. As the temperature rises, the pressure of the model rises above that 
of the free gas. We interpret this as a consequence of the decreasing mass function, Eq. 0, 
as one roughly expects that the average momentum of field excitations should increase with 
temperature. A difference between two kinds of modelling should be expected, since our model 
has positivity-violating quarks as its elementary constituents. 

Figure displays our results for the normalized trace anomaly at zero chemical potential. 
Notice that the model approaches the massless (conformally invariant) behavior at high 
temperature faster than the massive free gas. 
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The zero-temperature, finite chemical potential pressure (15) can be seen in Figure]^ Notice 
that the first excitations appear at a chemical potential p ~ Mthr = 0.467 GeV, a mass scale 
which is not explicitly present in the model action ([^. We thus interpret the threshold mass 
Mthr as a dynamically generated mass scale. Notice also that Mthr was used for the comparison 
between the model and a free gas in Figures and where a low-temperature agreement between 















Figure 1. Pressure for the model 
quarks (solid, red line), for bag model 
quarks (dashed, blue line, with B = 
(0.145 GeV)"^), and for a gas of free but 
massive quarks (dashed-dotted yellow line, 
with mass Mthr = 0.467 GeV). All lines are 
normalized by the Stefan-Boltzmann limit 
PSB = iTi'^NcNfT^/lSQ. 



Figure 2. . Trace anomaly at /x = 0 for 
dressed quarks (solid, red line), for massless 
quarks under a bag pressure (dashed, blue 
line, with B = (0.145 GeV)^), and for 
massive free quarks (dashed-dotted line, 
with mass = 0.467 GeV). 


the curves appears to be satisfactory. This indicates that a description purely in terms of a free 
gas of mass Mthr can only be reasonable just above the threshold, but not at higher densities. 
This is clearly a consequence of the momentum dependence of the mass function 0 - 

For completeness, we show in Figure]^ our results for the pressure at both finite temperature 
and finite chemical potential. As one could expect from general thermodynamical arguments, for 
a fixed value of chemical potential, increasing temperatures lead to higher pressures. Notice that, 
as soon as the temperature is nonzero, thermal activation allows field excitations for arbitrarily 
low chemical potentials. This can be seen from the nonzero values as well as a smoothening of 
the pressure curve already below ^ = Mthr- 

Let us finally notice that all our results are compatible with thermodynamical stability, a 
feature that is not always present in quark models with complex masses (for two studies, see, 

e.g., mm)- 

5. Summary 

In this work, we report on a first exploratory investigation of the thermodynamics of the confining 
quark model proposed in [T3j. We computed exactly the lowest-order partition function of the 
model at arbitrary temperature and chemical potential. Notice that, although the theory in its 
lowest nontrivial order is quadratic in the fields, it does not correspond to a free theory. Indeed, 
the quark propagator arising from the quadratic theory displays violation of reflection positivity 
and fits quite well lattice results for the quark mass function |24j . 

Our results for thermodynamical observables such as pressure and trace anomaly show that 
the microscopic confined degrees of freedom of the model do not correspond to a free gas of 
massive particles, even at lowest order. This can be understood from the dressing of quarks 
by the underlying gluons, which is encoded in the momentum-dependent mass function Q. In 
spite of the effective presence of particles with complex masses in the calculation of the partition 
function, we have not found any thermodynamical instabilities in the observables computed. 

























Figure 3. Zero-temperature limit of 
the pressure as a function of the chem¬ 
ical potential, normalized by the Stefan- 
Boltzmann limit Psb = /{12'k'^). 

Red, solid line: model result. Yellow, 
dotted-dashed line: free massive gas, with 
mass Mtfir = 0.467 GeV. 


Figure 4. Model pressure as a 

function of chemical potential for various 
temperatures. Solid, red: T = 0. Dotted, 
blue: T = 50 MeV. Dotted-dashed, yellow: 
T = 100 MeV. Solid with triangles, green: 
T = 150 MeV. Solid with circles, gray: 
T = 200 MeV. 
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